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LDA+DMFT is a novel computational technique for ab initio investigations of real materials 
with strongly correlated electrons, such as transition metals and their oxides. It combines the 
strength of conventional band structure theory in the local density approximation (LDA) with a 
modem many-body approach, the dynamical mean-field theory (DMFT). In the last few years 
LDA+DMFT has proved to be a powerful tool for the realistic modeling of strongly correlated 
electronic systems. In this paper the basic ideas and the set-up of the LDA+DMFT(X) approach, 
where X is the method used to solve the DMFT equations, are discussed. Results obtained with 
X=QMC (quantum Monte Carlo) and X=NCA (non-crossing approximation) are presented and 
compared. By means of the model system Lai-^Sr^TiOa we show that the method X matters 
qualitatively and quantitatively. Furthermore, we discuss recent results on the Mott-Hubbard 
metal-insulator transition in the transition metal oxide V2O3 and the a-7 transition in the 4f- 
electron system Ce. 
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1 Introduction 



The calculation of physical properties of electronic systems by controlled approximations 
is one of the most important challenges of modern theoretical solid state physics. In partic- 
ular, the physics of transition metal oxides - a singularly important group of materials both 
from the point of view of fundamental research and technological applications - may only 
be understood by explicit consideration of the strong effective interaction between the con- 
duction electrons in these systems. The investigation of electronic many-particle systems 
is made especially complicated by quantum statistics, and by the fact that the investigation 
of many phenomena require the application of non-perturbative theoretical techniques. 

From a microscopic point of view theoretical solid state physics is concerned with the 
investigation of interacting many-particle systems involving electrons and ions. However, 
it is an established fact that many electronic properties of matter are well described by the 
purely electronic Hamiltonian 



H = J d3r * + (r,cr) 
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where the crystal lattice enters only through an ionic potential. The applicability of this 
approach may be justified by the validity of the Born and Oppenheimer approximation. 1 
Here, ^ + (r, a) and \&(r, <r) are field operators that create and annihilate an electron at 
position r with spin a, A is the Laplace operator, m e the electron mass, e the electron 
charge, and 

VWrH- e 2 £r^7 V ee (r-r>) = £ £ * (2) 

i 1 11 r/r' 1 1 

denote the one-particle potential due to all ions i with charge eZi at given positions Rj, 
and the electron-electron interaction, respectively. 

While the ab initio Hamiltonian (1) is easy to write down it is impossible to solve 
exactly if more than a few electrons are involved. Numerical methods like Green's Func- 
tion Monte Carlo and related approaches have been used successfully for relatively modest 
numbers of electrons. Even so, however, the focus of the work has been on jellium and 
on light atoms and molecules like H, H2, 3 He, 4 He, see, e.g., the articles by Anderson, 
Bernu, Ceperley et al. in the present Proceedings of the NIC Winterschool 2002. Be- 
cause of this, one generally either needs to make substantial approximations to deal with 
the Hamiltonian (1), or replace it by a greatly simplified model Hamiltonian. At present 
these two different strategies for the investigation of the electronic properties of solids 
are applied by two largely separate groups: the density functional theory (DFT) and the 
many-body community. It is known for a long time already that DFT, together with its 
local density approximation (LDA), is a highly successful technique for the calculation 
of the electronic structure of many real materials. 2 However, for strongly correlated ma- 
terials, i.e., d- and /-electron systems which have a Coulomb interaction comparable to 
the band-width, DFT/LDA is seriously restricted in its accuracy and reliability. Here, the 



2 



model Hamiltonian approach is more general and powerful since there exist systematic 
theoretical techniques to investigate the many-electron problem with increasing accuracy. 
These many-body techniques allow one to describe qualitative tendencies and understand 
the basic mechanism of various physical phenomena. At the same time the model Hamil- 
tonian approach is seriously restricted in its ability to make quantitative predictions since 
the input parameters are not accurately known and hence need to be adjusted. One of the 
most successful techniques in this respect is the dynamical mean-field theory (DMFT) - 
a non-perturbative approach to strongly correlated electron systems which was developed 
during the past decade. 311 The LDA+DMFT approach, which was first formulated by 
Anisimov et al., 12 ' 13 combines the strength of DFT/LDA to describe the weakly correlated 
part of the ab initio Hamiltonian (1), i.e., electrons in s- and p-orbitals as well as the long- 
range interaction of the d- and /-electrons, with the power of DMFT to describe the strong 
correlations induced by the local Coulomb interaction of the d- or /-electrons. 

Starting from the ab initio Hamiltonian (1), the LDA+DMFT approach is presented in 
Section 2, including the DFT in Section 2.1, the LDA in Section 2.2, the construction of a 
model Hamiltonian in Section 2.3, and the DMFT in Section 2.4. As methods used to solve 
the DMFT we discuss the quantum Monte Carlo (QMC) algorithm in Section 2.5 and the 
non-crossing approximation (NCA) in Section 2.6. A simplified treatment for transition 
metal oxides is introduced in Section 2.7, and the scheme of a self-consistent LDA+DMFT 
in Section 2.8. As a particular example, the LDA+DMFT calculation for Lai-^Sr^TiOa 
is discussed in Section 3, emphasizing that the method X to solve the DMFT matters on 
a quantitative level. Our calculations for the Mott-Hubbard metal-insulator transition in 
V2O3 are presented in Section 4, in comparison to the experiment. Section 5 reviews our 
recent calculations of the Ce a-7 transition, in the perspective of the models referred to 
as Kondo volume collapse and Mott transition scenario. A discussion of the LDA+DMFT 
approach and its future prospects in Section 6 closes the presentation. 



2 The LDA+DMFT approach 
2.1 Density functional theory 

The fundamental theorem of DFT by Hohenberg and Kohn 14 (see, e.g., the review by 
Jones and Gunnarsson 2 ) states that the ground state energy is a functional of the elec- 
tron density which assumes its minimum at the ground state electron density. Following 
Levy, 15 this theorem is easily proved and the functional even constructed by taking the 
minimum (infimum) of the energy expectation value w.r.t. all (many-body) wave functions 
<p(ri<7i , ... tncttv) at a given electron number N which yield the electron density p(r): 

N 

E\p]=ini{(<p\H\<p) J (^<5(r-r0^>=p(r)}. (3) 

i=l 

However, this construction is of no practical value since it actually requires the eval- 
uation of the Hamiltonian (1). Only certain contributions like the Hartree energy 
^Hartree[p] = 3/ d 3 r' d 3 r V^r — r') p(y')p(y) and the energy of the ionic potential 
Eion [p] = j d 3 r Vion (r) p(r) can be expressed directly in terms of the electron density. 
This leads to 

E [p] = E kin [ P ] + E ion [p] + ^Hartrcc [p] + E xc [p] , (4) 
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where E^ [p] denotes the kinetic energy, and E xc [p] is the unknown exchange and cor- 
relation term which contains the energy of the electron-electron interaction except for the 
Hartree term. Hence all the difficulties of the many-body problem have been transferred 
into E xc [p] . While the kinetic energy £ki n cannot be expressed explicitly in terms of the 
electron density one can employ a trick to determine it. Instead of minimizing E[p] with 
respect to p one minimizes it w.r.t. a set of one-particle wave functions ipi related to p via 

N 

p(r)=^|^(r)| 2 . (5) 

i=i 

To guarantee the normalization of ipi, the Lagrange parameters are introduced such 
that the variation 5{E[p] + £.;[1 - / d 3 r\(pi(r)\ 2 }} /8<pi(r) = yields the Kohn-Sham 16 
equations: 

h2 a , tz ^ , fja j,r , SE *c\PY 



-A + Mon(r) + J d 3 r'V cc (r-r , )p(r') + 



■ ■ • , , • ... - <fii(r) = Si !fi(r). (6) 

2m e J op(r) J 

These equations have the same form as a one-particle Schrodinger equation which, a pos- 
teriori, justifies to calculate the kinetic energy by means of the one-particle wave-function 
ansatz. The kinetic energy of a one-particle ansatz which has the ground state density is, 
then, given by -EkinlPmin] = — J^iLii'fiil^ 2 ^ / (^ m e)\ l fii) if ^ e Pi are tne self-consistent 
(spin-degenerate) solutions of Eqs. (6) and (5) with lowest "energy" ej. Note, however, 
that the one-particle potential of Eq. (6), i.e., 

K ff (r) = Mon(r) + J d\'V cc (r-r')p(r') + , ( 7 ) 

is only an auxiliary potential which artificially arises in the approach to minimize E[p]. 
Thus, the wave functions <pi and the Lagrange parameters £j have no physical meaning at 
this point. Altogether, these equations allow for the DFT/LDA calculation, see the flow 
diagram Fig. 1. 



2.2 Local density approximation 

So far no approximations have been employed since the difficulty of the many-body prob- 
lem was only transferred to the unknown functional E xc [p] . For this term the local density 
approximation (LDA) which approximates the functional E xc [p] by a function that depends 
on the local density only, i.e., 

E xc [p}^ J d\ E^ A (p(r)), (8) 

was found to be unexpectedly successful. Here, E^ A (p(r)) is usually calculated from 
the perturbative solution 17 or the numerical simulation 18 of the jellium problem which is 
defined by Vi on (r) = const. 

In principle DFT/LDA only allows one to calculate static properties like the ground 
state energy or its derivatives. However, one of the major applications of LDA is the 
calculation of band structures. To this end, the Lagrange parameters Si are interpreted 
as the physical (one-particle) energies of the system under consideration. Since the true 
ground-state is not a simple one-particle wave-function, this is an approximation beyond 



4 



First principles information: 

atomic numbers, crystal structure (lattice, atomic positions) 



Choose initial electronic density p(r) 



Calculate effective potential using the LDA [Eq. (7)] 



Veff(r) = y ion (r) + / dV V ee (r - r')p(r') 



Sp(v) 



Solve Kohn-Sham equations [Eq. (6)] 

_|lv 2 + 14ff(r)- £i 
2m 



¥>i(r) = 



Calculate electronic density [Eq. (5)], 



JV 



p(r)^^|^(r)| 5 



Iterate to self-consistency 



Calculate band structure £i(k) [Eq. (6)], partial and total DOS, self-consistent 
Hamiltonian [Eq. (11)] LDA+DMFT, total energy E[p] [Eq. (3)], . . . 



Figure 1. Flow diagram of the DFT/LDA calculations. 



DFT. Actually, this approximation corresponds to the replacement of the Hamiltonian (1) 
by 



#lda = J2 j d3r 



+ 



2m 



A + V ion (r)+ j d 3 r' p(r')V cc (r-r') 

*(r,a). (9) 



^LDA 

Sp(r) 



For practical calculations one needs to expand the field operators w.r.t. a basis &u m , e.g., 
a linearized muffin-tin orbital (LMTO) 19 basis (i denotes lattice sites; I and m are orbital 
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indices). In this basis, 



* + (r,a)=^c^ m $ 4im (r) (10) 

ilm 

such that the Hamiltonian (9) reads 

H L DA = ^2 (Silm^l'm' Eilrn n^ lm +tu m jl'rn' C^ m Cj llm ,). (11) 



Here, n? lm — c^cf^, 



tamjJ'm' = (*« m - 7^ + ^ion(r) + j d 3 r'/,(r')F cc (r-r') + ^g^T 



Sp(r) 



®jl>m> 



(12) 



for iZm ^ j7'm' and zero otherwise; eu m denotes the corresponding diagonal part. 

As for static properties, the LDA approach based on the self-consistent solution of 
Hamiltonian (11) together with the calculation of the electronic density Eq. (5) [see the 
flow diagram Fig. 1] has also been highly successful for band structure calculations - 
but only for weakly correlated materials. 2 It is not reliable when applied to correlated 
materials and can even be completely wrong because it treats electronic correlations only 
very rudimentarily. For example, it predicts the antiferromagnetic insulator La2CuC>4 to be 
a non-magnetic metal 20 and also completely fails to account for the high effective masses 
observed in 4/-based heavy fermion compounds. 



2.3 Supplementing LDA with local Coulomb correlations 

Of prime importance for correlated materials are the local Coulomb interactions between 
d- and /-electrons on the same lattice site since these contributions are largest. This is due 
to the extensive overlap (w.r.t. the Coulomb interaction) of these localized orbitals which 
results in strong correlations. Moreover, the largest non-local contribution is the nearest- 
neighbor density-density interaction which, to leading order in the number of nearest- 
neighbor sites, yields only the Hartree term (see Ref. 4 and, also, Ref. 21) which is already 
taken into account in the LDA. To take the local Coulomb interactions into account, one 
can supplement the LDA Hamiltonian (11) with the local Coulomb matrix approximated 
by the (most important) matrix elements U™ m , (Coulomb repulsion and Z-component of 
Hund's rule coupling) and J mm ' (spin-flip terms of Hund's rule coupling) between the 
localized electrons (for which we assume i — id and I — 1^): 



H - i/LDA - #LDA + 2 X! Umm'^ilmcrnilm'cr' 

_1 \ ^ \ ^ j 4 ~t 

cy / j / j ^mm' c ilma C ilm'a C ilm'(7 C ilma' 



2 

i—i d .l—l d ma.m' 



(13) 



Here, the prime on the sum indicates that at least two of the indices of an operator have to be 
different, and a =\ (f) for a =f (J.). A term H^da * s subtracted to avoid double-counting 
of those contributions of the local Coulomb interaction already contained in Hlda- Since 
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there does not exist a direct microscopic or diagrammatic link between the model Hamil- 
tonian approach and LDA it is not possible to express H^ DA rigorously in terms of U, J 
and p. A commonly employed approximation for H^da assumes the LDA energy E^ DA 
ofH L %tobe 22 

^lda = ^Un d (n d - 1) - -J^n da {n ds - 1). (14) 

Here, n da = J2 m nu d rna = J2 m (hii d rna) is the total number of interacting electrons per 
spin, n d = J2a nd <*' U is the average Coulomb repulsion and J the average exchange or 
Hund's rule coupling. In typical applications we have Uj r j- m = U, J mm > = J, U™ m , = 
U — J — JSaa' for to 7^ to' (here, the first term J is due to the reduced Coulomb repulsion 
between different orbitals and the second term Jd aa > directly arises from the Z-component 
of Hund's rule coupling), and (with the number of interacing orbitals M) 

- _ U + (M - l)(U - J) + (M - 1)(U - 2J) 
2M - 1 

Since the one-electron LDA energies can be obtained from the derivatives of the total 
energy w.r.t. the occupation numbers of the corresponding states, the one-electron energy 
level for the non-interacting, paramagnetic states of (13) is obtained as 22 

Al dm = ^-^(^lda - E" DA ) = e lldm - U(n d - 1) + J -{n d - 1) (15) 

where Su d , n is defined in (11) and -Elda is the total energy calculated from Hlda (11). 
Furthermore we used n da = n d /2in the paramagnet. 

This leads to a new Hamiltonian describing the non-interacting system 

-^LDA = (fiilmJl'm'Alm™ilm+tilmJl'm'Cil m &jl' m ')> (^) 

ilmjl'm' ,a 

where e° ildTn is given by (15) for the interacting orbitals and e\ lm = eu, n for the non- 
interacting orbitals. While it is not clear at present how to systematically subtract H^ DA 
one should note that the subtraction of a Hartree-type energy does not substantially affect 
the overall behavior of a strongly correlated paramagnetic metal in the vicinity of a Mott- 
Hubbard metal-insulator transition (see also Section 2.7). 

In the following, it is convenient to work in reciprocal space where the matrix elements 
of H^ DA , i.e., the LDA one-particle energies without the local Coulomb interaction, are 
given by 



(-^LDAWVm^'Z'm' — (HhDA (k) )qlm,q> I'm' 



Vqlm,q'l'm'Oql,q d l d 



U{n d -\)- J -{n d -l) 



(17) 



Here, q is an index of the atom in the elementary unit cell, (Hu)A(^))qim,q'l'm' is the 
matrix element of (11) in k-space, and q d denotes the atoms with interacting orbitals in 
the unit cell. The non-interacting part, #l DA , supplemented with the local Coulomb in- 
teraction forms the (approximated) ab initio Hamiltonian for a particular material under 
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investigation: 

i—i^A—ld mer,m' 'a' 

-- V V' 7 ,c f c f c c (18) 

2 / j / j mm'^ilma ilm'cr ilm'a ilmcr \ iv -v 
i—id^—ld vaa.ra' 

To make use of this ab initio Hamiltonian it is still necessary to determine the Coulomb 
interaction U. To this end, one can calculate the LDA ground state energy for different 
numbers of interacting electrons rid ("constrained LDA" 23 ) and employ Eq. (14) whose 
second derivative w.r.t. rid yields U. However, one should keep in mind that, while the 
total LDA spectrum is rather insensitive to the choice of the basis, the calculation of U 
strongly depends on the shape of the orbitals which are considered to be interacting. E.g., 
for LaTi0 3 at a Wigner Seitz radius of 2.37 a.u. for Ti a LMTO-ASA calculation 24 using 
the TB-LMTO-ASA code 19 yielded U = 4.2 eV in comparison to the value U = 3.2 eV 
calculated by ASA-LMTO within orthogonal representation. 25 Thus, an appropriate basis 
like LMTO is mandatory and, even so, a significant uncertainty in U remains. 



2.4 Dynamical mean-field theory 

The many-body extension of LDA, Eq. (18), was proposed by Anisimov et al. 22 in the 
context of their LDA+U approach. Within LDA+U the Coulomb interactions of (18) are 
treated within the Hartree-Fock approximation. Hence, LDA+U does not contain true 
many-body physics. While this approach is successful in describing long-range ordered, 
insulating states of correlated electronic systems it fails to describe strongly correlated 
paramagnetic states. To go beyond LDA+U and capture the many-body nature of the 
electron-electron interaction, i.e., the frequency dependence of the self-energy, various ap- 
proximation schemes have been proposed and applied recently. 12 ' 26-30 One of the most 
promising approaches, first implemented by Anisimov et al., 12 is to solve (18) within 
DMFT 311 ("LDA+DMFT"). Of all extensions of LDA only the LDA+DMFT approach 
is presently able to describe the physics of strongly correlated, paramagnetic metals with 
well-developed upper and lower Hubbard bands and a narrow quasiparticle peak at the 
Fermi level. This characteristic three-peak structure is a signature of the importance of 
many-body effects. 7 ' 8 

During the last ten years, DMFT has proved to be a successful approach to investigate 
strongly correlated systems with local Coulomb interactions. 11 It becomes exact in the 
limit of high lattice coordination numbers 3 ' 4 and preserves the dynamics of local interac- 
tions. Hence, it represents a dynamical mean-field approximation. In this non-perturbative 
approach the lattice problem is mapped onto an effective single-site problem (see Fig. 2) 
which has to be determined self-consistently together with the k-integrated Dyson equation 
connecting the self energy X and the Green function G at frequency u>: 

G qlm , q . Vm .{w)=-2 r /rf 3 fc([^l + M l-i70 DA (k)-SH]" 1 ) . (19) 

\>B J ^ ' qlm,q'l'm' 

Here, 1 is the unit matrix, /j, the chemical potential, the matrix i7° DA (k) is defined in 
(17), denotes the self-energy matrix which is non-zero only between the interacting 
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Figure 2. If the number of neighboring lattice sites goes to infinity, the central limit theorem holds and fluctua- 
tions from site-to-site can be neglected. This means that the influence of these neighboring sites can be replaced 
by a mean influence, the dynamical mean-field described by the self energy Y,f m (u>). This DMFT problem is 
equivalent to the self-consistent solution of the k-integrated Dyson equation (21) and the multi-band Anderson 
impurity model Eq. (20). 

orbitals, implies the inversion of the matrix with elements n (=qlm), n'(=q'l'm'), 

and the integration extends over the Brillouin zone with volume Vb ■ 

The DMFT single-site problem depends on Q^)^ 1 = G(w) _1 + S(w) and is 
equivalent 7 ' 8 to an Anderson impurity model (the history and the physics of this model 
is summarized by Anderson in Ref. 31) if its hybridization A(w) satisfies = 
u> — J du' A(co') / (u> — uj 1 ). The local one-particle Green function at a Matsubara fre- 
quency iu v = i(2v + (/?: inverse temperature), orbital index m (I = l d , q = q d ), 
and spin a is given by the following functional integral over Grassmann variables ip and 

Here, Z = J P[i/)]D[f ]f m C exp(.4[V>, ip* , G' 1 ]) is the partition function and the 
single-site action A has the form (the interaction part of A is in terms of the "imaginary 
time" r, i.e., the Fourier transform of w„) 

A&,r,g- 1 } = E vc&rv^ 

is,<7.m 

3 

~E' U ™' I dT ^{rW m (r)C>(r)C' (r) 

m<j,ma' q 

S3 

+ \ E' J mm'[ drr m *(rW m (rW m *(rW^(r) . (21) 

raa.ra q 

This single-site problem (20) has to be solved self-consistently together with the It- 
integrated Dyson equation (19) to obtain the DMFT solution of a given problem, see the 
flow diagram Fig. 3. 

Due to the equivalence of the DMFT single-site problem and the Anderson impurity 
problem a variety of approximative techniques have been employed to solve the DMFT 
equations, such as the iterated perturbation theory (IPT) 7 ' 1 1 and the non-crossing approxi- 
mation (NCA), 32-34 as well as numerical techniques like quantum Monte Carlo simulations 
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Choose an initial self-energy £ 




Calculate G from E via the k-integrated Dyson Eq. (19): 

G q i m , q , l >m>(uj)=^ r f ^([c.l + Ml-I/VW-SMr 1 ) , ,„ , 

VB J \ ' qlm,q'l'm' 




g = (G- 1 + ^)- 1 




Calculate G from Q via the DMFT single-site problem Eq. (20) 




Snow = G^ 1 - G^ 1 


Iterate with S = E now until convergence, i.e. S — S nC w < e 



Figure 3. Flow diagram of the DMFT self-consistency cycle. 



(QMC), 35 exact diagonalization (ED), 11 - 36 or numerical renormalization group (NRG). 37 
QMC and NCA will be discussed in more detail in Section 2.5 and 2.6, respectively. IPT is 
non-self-consistent second-order perturbation theory in U for the Anderson impurity prob- 
lem (20) at half-filling. It represents an ansatz that also yields the correct perturbational 
£/ 2 -term and the correct atomic limit for the self-energy off half-filling, 38 for further details 
see Refs. 12,26, 38. ED directly diagonalizes the Anderson impurity problem at a limited 
number of lattice sites and orbitals. NRG first replaces the conduction band by a discrete 
set of states at DA - ™ (D: bandwidth; n = 0, Af s ) and then diagonalizes this problem 
iteratively with increasing accuracy at low energies, i.e., with increasing Af s . In principle, 
QMC and ED are exact methods, but they require an extrapolation, i.e., the discretization 
of the imaginary time At — > (QMC) or the number of lattice sites of the respective 
impurity model N s — > oo (ED), respectively. 

In the context of LDA+DMFT we refer to the computational schemes to solve the 
DMFT equations discussed above as LDA+DMFT(X) where X=IPT, 12 NCA, 30 QMC 24 
have been investigated in the case of Lai-^Sr^TiOa. The same strategy was formulated 
by Lichtenstein and Katsnelson 26 as one of their LDA++ approaches. Lichtenstein and 
Katsnelson applied LDA+DMFT(IPT), 42 and were the first to use LDA+DMFT(QMC), 43 
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to investigate the spectral properties of iron. Recently, also V2O3, Ca2- x Sr x Ru04, ' 6 
Ni, 47 Fe, 47 Pu, 48 - 49 and Ce 50 - 51 have been studied by LDA+DMFT. Realistic investiga- 
tions of itinerant ferromagnets (e.g., Ni) have also recently become possible by combining 
density functional theory with multi-band Gutzwiller wave functions. 52 



2.5 QMC method to solve DMFT 

The self-consistency cycle of the DMFT (Fig. 3) requires a method to solve for the dynam- 
ics of the single-site problem of DMFT, i.e., Eq. (20). The QMC algorithm by Hirsch and 
Fye 35 is a well established method to find a numerically exact solution for the Anderson 
impurity model and allows one to calculate the impurity Green function G at a given Q~ x 
as well as correlation functions. In essence, the QMC technique maps the interacting elec- 
tron problem Eq. (20) onto a sum of non-interacting problems where the single particle 
moves in a fluctuating, time-dependent field and evaluates this sum by Monte Carlo sam- 
pling, see the flow diagram Fig. 4 for an overview. To this end, the imaginary time interval 
[0,/3] of the functional integral Eq. (20) is discretized into A steps of size At = /3/A, 
yielding support points 77 = I At with I = 1 ... A. Using this Trotter discretization, the 
integral dr is transformed to the sum J2t=i Ar tne ex P onen tial terms in Eq. (20) 
can be separated via the Trotter-Suzuki formula for operators A and B 5i 

A 

e -0(A+B) = JJ e -Arl e -ArS + 0{At) ^ (22) 
2 = 1 

which is exact in the limit At — > 0. The single site action A of Eq. (21) can now be written 
in the discrete, imaginary time as 

A-l 

A[^,r,G- 1 ]=Ar 2 ]T r m i*G a m -\lAr-l'ArW m v 

aml,l'=0 

, A-l 

"9 Ar E , (23) 

1=0 

where the first term was Fourier-transformed from Matsubara frequencies to imaginary 
time. In a second step, the M(2M — 1) interaction terms in the single site action A are 
decoupled by introducing a classical auxiliary field sf£ m , : 

exp {^U™ m ,(r ml *r m i - r m nV m nf 

\ J2 ^p{ A ^rw^w(C ; >™/-^r^/)}, (24) 

a""' , =±1 

where cosh( Af^ m , ) = exp( AtU^, /2) and M is the number of interacting orbitals. This 
so-called discrete Hirsch-Fye-Hubbard-Stratonovich transformation can be applied to the 
Coulomb repulsion as well as the Z-component of Hund's rule coupling. 54 It replaces the 
interacting system by a sum of AM (2M — 1) auxiliary fields sf£ m , . The functional integral 
can now be solved by a simple Gauss integration because the Fermion operators only enter 
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quadratically, i.e., for a given configuration s = {s^ lm ,} of the auxiliary fields the system 
is non-interacting. The quantum mechanical problem is then reduced to a matrix problem 

g| W2 = ^E E E [( M i s r\ l2 U detM ^ (25) 

I m' a' .m" a" s a " f' , — ±1 m<T 

with the partition function Z, the matrix 

M|f = Ar 2 ^- 1 + Je"*" + 1 - e~" x " (26) 
and the elements of the matrix A£f 

\crs c \ \crcr' ^cc' £ . crcr ' ^OT* 

A mJi' — °«' /-^i A mm' a mm' s lmm'- \ L 1 > 

m'a' 

Here o™ m i = 2Q(a' — a + S aa >[m' — m] — 1) changes sign if (mcr) and (m'a 1 ) are 
exchanged. For more details, e.g., for a derivation of Eq. (26) for the matrix M, see 
Refs. 11,35. 

Since the sum in Eq. (25) consists of 2 AM ( 2M_1 ) addends, a complete summation 
for large A is computationally impossible. Therefore the Monte Carlo method, which is 
often an efficient way to calculate high-dimensional sums and integrals, is employed for 
importance sampling of Eq. (25). In this method, the integrand F(x) is split up into a 
normalized probability distribution P and the remaining term O: 

J dxF(x) = J dxO(x)P(x) = (0) P (28) 

with 

dxP(x) = 1 and P{x) > 0. (29) 



In statistical physics, the Boltzmann distribution is often a good choice for the function P: 

P(x) = jexp(-l3E(x)). (30) 
For the sum of Eq. (25), this probability distribution translates to 

p ( s )=^ri dctM - < 3i > 



z 

ma 



with the remaining term 



0(sf mlll2 = [(Mg)-\ h . (32) 

Instead of summing over all possible configurations, the Monte Carlo simulation gen- 
erates configurations Xj with respect to the probability distribution P(x) and averages the 
observable O(x) over these Xj. Therefore the relevant parts of the phase space with a large 
Boltzmann weight are taken into account to a greater extent than the ones with a small 
weight, coining the name importance sampling for this method. With the central limit 
theorem one gets for J\f statistically independent addends the estimate 

^ p = Jf E 0{x l )±-^=^j(0*) P -(0) P . (33) 
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Choose random auxiliary field configuration s = {s^ m ,} 



Calculate the current Green function G cur from Eq. (32) 

(Gcur)'rhl 1 l 2 = [(M^n) 1 ]; 1 ; 2 

with M fromEq. (26) and the input G^v)' 1 = G a m {u v )- X + T, a m {w v ). 



Do NWU times (warm up sweeps) 



MC-sweep (G cur ,s) 



Do NMC times (measurement sweeps) 



MC-sweep (G cur ,s) 



G = G + G cur /NMC 



Figure 4. Flow diagram of the QMC algorithm to calculate the Green function matrix G using the procedure 
MC-sweep of Fig. 5. 



Choose M(2M - 1) A times a set (I mm' a <r'), 

define s now to be s except for the element sf r f lm , which has opposite sign. 




Calculate flip probability P s ^ Snow = min{l, P(s 


aew)AP(s)} With 




P(s ncw )/P(s) = J] dot / [] det M£ 

m.rr mrr 




and M from Eq. (26). 






Random number e (0, 1) < P s ^s„ ow ? ^^^^ 




yes 


^^^^ no 




s = s ncw ; recalculate G cur according to Eq. (32). 


Keep s 



Figure 5. Procedure MC-sweep using the Metropolis 55 rule to change the sign of s^ m , . The recalculation of 
Gcur, i.e., the matrix M of Eq. (26), simplifies to C(A 2 ) operations if only one s" a , changes sign. 11,35 
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Here, the error and with it the number of needed addends J\f is nearly independent of 
the dimension of the integral. The computational effort for the Monte Carlo method is 
therefore only rising polynomially with the dimension of the integral and not exponentially 
as in a normal integration. Using a Markov process and single spin-flips in the auxiliary 
fields, the computational cost of the algorithm in leading order of A is 

2aM(2M - 1)A 3 x number of MC-sweeps, (34) 

where a is the acceptance rate for a single spin-flip. 

The advantage of the QMC method (for the algorithm see the flow diagram Fig. 4) is 
that it is (numerically) exact. It allows one to calculate the one-particle Green function as 
well as two-particle (or higher) Green functions. On present workstations the QMC ap- 
proach is able to deal with up to seven interacting orbitals and temperatures above about 
room temperature. Very low temperatures are not accessible because the numerical effort 
grows like A 3 oc 1/T 3 . Since the QMC approach calculates G(t) or G(iu n ) with a sta- 
tistical error, it also requires the maximum entropy method 56 to obtain the Green function 
G(lu) at real (physical) frequencies lo. 



2.6 NCA method to solve DMFT 

The NCA approach is a resolvent perturbation theory in the hybridization parameter A(w) 
of the effective Anderson impurity problem. 32 Thus, it is reliable if the Coulomb interac- 
tion U is large compared to the band-width and also offers a computationally inexpensive 
approach to check the general spectral features in other situations. 

To see how the NCA can be adapted for the DMFT, let us rewrite Eq. (19) as 

= w k £ [z " H ° LDA(k) 1 (35) 

k 

where z = u) + iQ + + [i. Again, Hl DA (k), and hence G"(C) and G„(z) are matrices 
in orbital space. Note that has nonzero entries for the correlated orbitals only. 
On quite general grounds, Eq. (35) can be cast into the form 



where 



E ° = ^- k Y, H °LDA(V (37) 

with the number of k points Nk and 

lim $te{A a (uj + i6)} = . (38) 

Given the the matrix E°, the Coulomb matrix U and the hybridization matrix A a (z), 
we are now in a position to set up a resolvent perturbation theory with respect to A a (z). 
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To this end, we first have to diagonalize the local Hamiltonian 

H local = ^""^ ^""^ c \lma^qlm,q'l'm' c qlma 
cr qml q'm'V 



^2 ^mm' n qdldma n q d l d m' a' 

mcr m' cr' 

(39) 

2 / ^ / j rnm qdldma qdldm'a qdldm'a qdldma 



ma m' 



= ^E a \a)(a\ 



with local eigenstates \a) and energies E a . In contrast to the QMC, this approach allows 
one to take into account the full Coulomb matrix plus spin-orbit coupling. 

With the states \a) defined above, the fermionic operators with quantum numbers n = 
(q, I, m) are expressed as 



cL =£(£>&)»</ji, 



(40) 

a, 

The key quantity for the resolvent perturbation theory is the resolvent R(z), which obeys 
a Dyson equation 32 

R(z) = R°(z) + R°(z)S(z)R(z) , (41) 

where i?° /3 (z) = l/(z — E a )5 a p and S a p(z) denotes the self-energy for the local states 
due to the coupling to the environment through A(z). 

Retaining only the lowest-, i.e. 2" d -order terms leads to a set of self-consistent integral 
equations 



The self-energy S a p{z) can be expressed as power series in the hybridization A(z). 32 



s m = EEE / v /(£) ( DK JSK K 'm a 'fS'(z+e)D^ 

a kk ! a.' (3' 

+ EEE/ f (WOO) D% a r?'{e)R a , l) ,{z-e) (d^J 

a kk' a.' 0' 



(42) 



to determine S a p(z), where /(e) denotes Fermi's function and T(e) = 
— 3m {A(e + i0 + )}. The set of equations (42) are in the literature referred to as 
non-crossing approximation (NCA), because, when viewed in terms of diagrams, they 
contain no crossing of band-electron lines. In order to close the cycle for the DMFT, we 
still have to calculate the true local Green function G a (z). This, however, can be done 
within the same approximation with the result 

Gr'M-^EE^-'(^-')Y^^^ (z) ^ v ' (z+H ■ (43) 



a, a' v.v' 
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Here, Z\ oca \ = \ <X> — - — : — R aa (z) denotes the local partition function and is the 

a J 

inverse temperature. 

Like any other technique, the NCA has its merits and disadvantages. As a self- 
consistent resummation of diagrams it constitutes a conserving approximation to the An- 
derson impurity model. Furthermore, it is a (computationally) fast method to obtain dy- 
namical results for this model and thus also within DMFT. However, the NCA is known to 
violate Fermi liquid properties at temperatures much lower than the smallest energy scale 
of the problem and whenever charge excitations become dominant. 34 ' 57 Hence, in some 
parameter ranges it fails in the most dramatic way and must therefore be applied with 
considerable care. 34 



2.7 Simplifications for transition metal oxides with well separated e g - and 
t 2g -bands 

Many transition metal oxides are cubic perovskites, with only a slight distortion of the 
cubic crystal structure. In these systems the transition metal c?-orbitals lead to strong 
Coulomb interactions between the electrons. The cubic crystal-field of the oxygen causes 
the d-orbitals to split into three degenerate t 2g - and two degenerate e g -orbitals. This split- 
ting is often so strong that the t 2g - or e g -bands at the Fermi energy are rather well sepa- 
rated from all other bands. In this situation the low-energy physics is well described by 
taking only the degenerate bands at the Fermi energy into account. Without symmetry 
breaking, the Green function and the self-energy of these bands remain degenerate, i.e., 

G q l m . q il'm'{z) = G(z)8 q l m . q i p m > and 51qlm,q'l'm'(z) = ^(z)8qlm,q'l'm' f° r I = Id 

q = <id (where Id and qa denote the electrons in the interacting band at the Fermi energy). 
Downfolding to a basis with these degenerate ^-Z^-bands results in an effective Hamilto- 
nian (where indices I — Id and g = ^ are suppressed) 

G mm ,(Lo) = ^-Jd 3 k (M + Ml-^K^-SH]" 1 )^,. (44) 

Due to the diagonal structure of the self-energy the degenerate interacting Green function 
can be expressed via the non-interacting Green function G°(uj): 

GM-GV-SM) = / de N ° {€ ] . (45) 

J w-S(w)-e 

Thus, it is possible to use the Hilbert transformation of the unperturbed LDA-calculated 
density of states (DOS) N°(e), i.e., Eq. (45), instead of Eq. (19). This simplifies the 
calculations considerably. With Eq. (45) also some conceptual simplifications arise: (i) 
the subtraction of -Hlda m (45) only results in an (unimportant) shift of the chemical 
potential and, thus, the exact form of H^da i s irrelevant; (ii) Luttinger's theorem of Fermi 
pinning holds, i.e., the interacting DOS at the Fermi energy is fixed at the value of the non- 
interacting DOS at T = within a Fermi liquid; (iii) as the number of electrons within the 
different bands is fixed, the LDA+DMFT approach is automatically self-consistent. 

In this context it should be noted that the approximation Eq. (45) is justified only if the 
overlap between the t 2g orbitals and the other orbitals is rather weak. 
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2.8 Extensions of the LDA+DMFT scheme 



In the present form of the LDA+DMFT scheme the band-structure input due to LDA and 
the inclusion of the electronic correlations by DMFT are performed as successive steps 
without subsequent feedback. In general, the DMFT solution will result in a change of the 
occupation of the different bands involved. This changes the electron density p(r) and, 
thus, results in a new LDA-Hamiltonian Hlda (H) since Hlda depends on p(r). At 
the same time also the Coulomb interaction U changes and needs to be determined by a 
new constrained LDA calculation. In a self-consistent LDA+DMFT scheme, #lda and U 
would define a new Hamiltonian (18) which again needs to be solved within DMFT, etc., 
until convergence is reached: 

DMFT 

PO) - #lda, U n Um p(r) 

t I (46) 



Without Coulomb interaction (U = 0) this scheme reduces to the self-consistent solution of 
the Kohn-Sham equations. A self-consistency scheme similar to Eq. (46) was employed by 
Savrasov and Kotliar 49 in their calculation of Pu. An ab initio DMFT scheme formulated 
directly in the continuum was recently proposed by Chitra and Kotliar. 58 



3 Comparison of different methods to solve DMFT: the model system 
Lai^Sr^TiOs 

The stoichiometric compound LaTiC>3 is a cubic perovskite with a small orthorhombic 
distortion (Z Ti — O — Ti w 155 ) 59 and is an antiferromagnetic insulator 60 below 
Tn = 125 K. 61 Above TV, or at low Sr-doping x, and neglecting the small orthorhombic 
distortion (i.e., considering a cubic structure with the same volume), LaTiOa is a strongly 
correlated, but otherwise simple paramagnet with only one 3d-electron on the trivalent Ti 
sites. This makes the system a perfect trial candidate for the LDA+DMFT approach. 

The LDA band-structure calculation for undoped (cubic) LaTiC>3 yields the DOS 
shown in Fig. 6 which is typical for early transition metals. The oxygen bands, rang- 
ing from — 8.2 eV to — 4.0 eV, are filled such that Ti is three- valent. Due to the crystal-field 
splitting, the Ti 3d-bands separates into two empty e s -bands and three degenerate ti g - 
bands. Since the i 2ff -bands at the Fermi energy are well separated also from the other bands 
we employ the approximation introduced in section 2.5 which allows us to work with the 
LDA DOS [Eq. (45)] instead of the full one-particle Hamiltonian #£ DA of [Eq. (19)]. In 
the LDA+DMFT calculation, Sr-doping x is taken into account by adjusting the chemical 
potential to yield n = 1 — x = 0.94 electrons within the i 2g -bands, neglecting effects 
disorder and the x-dependence of the LDA DOS (note, that Sr and Ti have a very similar 
band structure within LDA). There is some uncertainty in the LDA-calculated Coulomb 
interaction parameter U <~ 4 — 5 eV (for a discussion see Ref. 24) which is here assumed 
to be spin- and orbital-independent. In Fig. 7, results for the spectrum of Lao.94Sro.o6Ti03 
as calculated by LDA+DMFT(IPT, NCA, QMC) for the same LDA DOS at T w 1000 K 
and U = 4 eV are compared. 24 In Ref. 24 the formerly presented IPT 12 and NCA 30 spec- 
tra were recalculated to allow for a comparison at exactly the same parameters. All three 
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Energy, eV 

Figure 6. Densities of states of LaTiC>3 calculated with LDA-LMTO. Upper figure: total DOS; lower figure: 
partial ti g (solid lines) and e g (dashed lines) DOS [reproduced from Ref.24]. 



methods yield the typical features of strongly correlated metallic paramagnets: a lower 
Hubbard band, a quasi-particle peak (note that IPT produces a quasi-particle peak only 
below about 250K which is therefore not seen here), and an upper Hubbard band. By 
contrast, within LDA the correlation-induced Hubbard bands are missing and only a broad 
central quasi-particle band (actually a one-particle peak) is obtained (Fig. 6). 

While the results of the three evaluation techniques of the DMFT equations (the ap- 
proximations IPT, NCA and the numerically exact method QMC) agree on a qualitative 
level, Fig. 7 reveals considerable quantitative differences. In particular, the IPT quasi- 
particle peak found at low temperatures (see right inset of Fig. 7) is too narrow such that 
it disappears already at about 250 K and is, thus, not present at T w 1000 K. A similarly 
narrow IPT quasi-particle peak was found in a three-band model study with Bethe-DOS 
by Kajueter and Kotliar. 38 Besides underestimating the Kondo temperature, IPT also pro- 
duces notable deviations in the shape of the upper Hubbard band. Although NCA comes 
off much better than IPT it still underestimates the width of the quasiparticle peak by a 
factor of two. Furthermore, the position of the quasi-particle peak is too close to the lower 
Hubbard band. In the left inset of Fig. 7, the spectra at the Fermi level are shown. At the 
Fermi level, where at sufficiently low temperatures the interacting DOS should be pinned 
at the non-interacting value, the NCA yields a spectral function which is almost by a fac- 
tor of two too small. The shortcomings of the NCA-results, with a too small low-energy 
scale and too much broadened Hubbard bands for multi-band systems, are well understood 
and related to the neglect of exchange type diagrams. 63 Similarly, the deficiencies of the 
IPT-results are not entirely surprising in view of the semi-phenomenological nature of this 
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Figure 7. Spectrum of Lao.94Sr .o6Ti0 3 as calculated by LDA+DMFT(X) at T = 0.1 eV (ss 1000 K) and 
U = 4 eV employing the approximations X=IPT, NCA, and numerically exact QMC. Inset left: Behavior at 
the Fermi level including the LDA DOS. Inset right: X=IPT and NCA spectra at T = 80 K [reproduced from 
Ref.24]. 



approximation, especially for a system off half filling. 

This comparison shows that the choice of the method used to solve the DMFT equations 
is indeed important, and that, at least for the present system, the approximations 1PT and 
NCA differ quantitatively from the numerically exact QMC. Nevertheless, the NCA gives 
a rather good account of the qualitative spectral features and, because it is fast and can 
often be applied to comparatively low temperatures, can yield an overview of the physics 
to be expected. 

Photoemission spectra provide a direct experimental tool to study the electronic struc- 
ture and spectral properties of electronically correlated materials. A comparison of 
LDA+DMFT(QMC) at 1000 K 65 with the experimental photoemission spectrum 64 of 
Lao.94Sro.o6Ti03 is presented in Fig 8. To take into account the uncertainty in U, 24 we 
present results for U = 3.2, 4.25 and 5 eV. All spectra are multiplied with the Fermi step 
function and are Gauss-broadened with a broadening parameter of 0.3 eV to simulate the 
experimental resolution. 64 LDA band structure calculations, the results of which are also 
presented in Fig. 8, clearly fail to reproduce the broad band observed in the experiment at 1- 
2 eV below the Fermi energy. 64 Taking the correlations between the electrons into account, 
this lower band is easily identified as the lower Hubbard band whose spectral weight orig- 
inates from the quasi-particle band at the Fermi energy and which increases with U. The 
best agreement with experiment concerning the relative intensities of the Hubbard band 
and the quasi-particle peak and, also, the position of the Hubbard band is found for U = 5 
eV. The value U — 5 eV is still compatible with the ab initio calculation of this parameter 
within LDA. 24 One should also bear in mind that photoemission experiments are sensitive 
to surface properties. Due to the reduced coordination number at the surface the bandwidth 
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Energy (eV) 

Figure 8. Comparison of the experimental photoemission spectrum, 64 the LDA result, and the 
LDA+DMFT(QMC) calculation for Lao.94Sro.o6TiC>3 (i.e., 6% hole doping) and different Coulomb interaction 
U = 3.2, 4.25, and 5 eV [reproduced from Ref.24]. 

is likely to be smaller, and the Coulomb interaction less screened, i.e., larger. Both effects 
make the system more correlated and, thus, might also explain why better agreement is 
found for U = 5 eV. Besides that, also the polycrystalline nature of the sample, as well 
as spin and orbital 66 fluctuation not taken into account in the LDA+DMFT approach, will 
lead to a further reduction of the quasi-particle weight. 

4 Mott-Hubbard metal-insulator transition in V 2 3 

One of the most famous examples of a cooperative electronic phenomenon occurring 
at intermediate coupling strengths is the transition between a paramagnetic metal and a 
paramagnetic insulator induced by the Coulomb interaction between the electrons - the 
Mott-Hubbard metal-insulator transition. The question concerning the nature of this tran- 
sition poses one of the fundamental theoretical problems in condensed matter physics. 67 
Correlation-induced metal-insulator transitions (MIT) are found, for example, in transition 
metal oxides with partially filled bands near the Fermi level. For such systems band theory 
typically predicts metallic behavior. The most famous example is V2O3 doped with Cr as 
shown in Fig. 9. While at low temperatures V2O3 is an antiferromagnetic insulator with 
monoclinic crystal symmetry, it has a corundum structure in the high-temperature param- 
agnetic phase. All transitions shown in the phase diagram are of first order. In the case 
of the transitions from the high-temperature paramagnetic phases into the low-temperature 
antiferromagnetic phase this is naturally explained by the fact that the transition is accom- 
panied by a change in crystal symmetry. By contrast, the crystal symmetry across the MIT 
in the paramagnetic phase remains intact, since only the ratio of the c/ a axes changes dis- 
continuously. This may be taken as an indication for the predominantly electronic origin 
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Figure 9. Experimental phase diagram of V2O3 doped with Cr and Ti [reproduced from Ref. 68]. Doping V2O3 
effects the lattice constants in a similar way as applying pressure (generated either by a hydrostatic pressure P, or 
by changing the V-concentration from V2O3 to V2- M 03) and leads to a Mott-Hubbard transition between the 
paramagnetic insulator (PI) and metal (PM). At lower temperatures, a Mott-Heisenberg transition between the 
paramagnetic metal (PM) and the antiferromagnetic insulator (AFI) is observed. 



of this transition which is not accompanied by any conventional long-range order. From a 
models point of view the MIT is triggered by a change of the ratio of the Coulomb interac- 
tion U relative to the bandwidth W . Originally, Mott considered the extreme limits W = 
(when atoms are isolated and insulating) and {7 = where the system is metallic. While it 
is simple to describe these limits, the crossover between them, i.e., the metal-insulator tran- 
sition itself, poses a very complicated electronic correlation problem. Among others, this 
metal-insulator transition has been addressed by Hubbard in various approximations 69 and 
by Brinkman and Rice within the Gutzwiller approximation. 70 During the last few years, 
our understanding of the MIT in the one-band Hubbard model has considerably improved, 
in particular due to the application of dynamical mean-field theory. 71 

Both the paramagnetic metal V2O3 and the paramagnetic insulator 
(V0.962CT0. 038)203 have the same corundum crystal structure with only slightly 
different lattice parameters. 72,73 Nevertheless, within LDA both phases are found to be 
metallic (see Fig. 10). The LDA DOS shows a splitting of the five Vanadium d-orbitals 
into three t2 9 states near the Fermi energy and two e a g states at higher energies. This 
reflects the (approximate) octahedral arrangement of oxygen around the vanadium atoms. 
Due to the trigonal symmetry of the corundum structure the t2 g states are further split 
into one a\ g band and two degenerate e g bands, see Fig. 10. The only visible difference 
between (Vo.962Cr . 038)203 and V2O3 is a slight narrowing of the t 2g and e a g bands by 
s» 0.2 and 0.1 eV, respectively as well as a weak downshift of the centers of gravity of 
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Figure 10. Left: Scheme of the 3d levels in the corundum crystal structure. Right: Partial LDA DOS of the 3d 
bands for paramagnetic metallic V2O3 and insulating (Vo. 962 Cro. 038)203 [reproduced from Ref.44]. 

both groups of bands for V 2 3 . In particular, the insulating gap of the Cr-doped system is 
seen to be missing in the LDA DOS. Here we will employ LDA+DMFT(QMC) to show 
explicitly that the insulating gap is caused by the electronic correlations. In particular, we 
make use of the simplification for transition metal oxides described in Section 2.7 and 
restrict the LDA+DMFT(QMC) calculation to the three ti g bands at the Fermi energy, 
separated from the e a g and oxygen bands. 

While the Hund's rule coupling J is insensitive to screening effects and may, thus, 
be obtained within LDA to a good accuracy (J = 0.93 eV 25 ), the LDA-calculated value 
of the Coulomb repulsion U has a typical uncertainty of at least 0.5 eV. 24 To overcome 
this uncertainty, we study the spectra obtained by LDA+DMFT(QMC) for three different 
values of the Hubbard interaction (U = 4.5, 5.0, 5.5) in Fig. 11. All QMC results presented 
were obtained for T = 0.1 eV. However, simulations for V2O3 at U = 5 eV, T = 0.143 
eV, and T = 0.067 eV suggest only a minor smoothing of the spectrum with increasing 
temperature. From the results obtained we conclude that the critical value of U for the 
MIT is at about 5 eV: At U = 4.5 eV one observes pronounced quasiparticle peaks at 
the Fermi energy, i.e., characteristic metallic behavior, even for the crystal structure of the 
insulator (Vo.962Cr . 038)203, while at U = 5.5 eV the form of the calculated spectral 
function is typical for an insulator for both sets of crystal structure parameters. At U = 5.0 
eV one is then at, or very close to, the MIT since there is a pronounced dip in the DOS at the 
Fermi energy for both a\ g and e£ orbitals for the crystal structure of (Vo.962Cr . 038)203, 
while for pure V2O3 one still finds quasiparticle peaks. (We note that at T w 0.1 eV one 
only observes metallic-like and insulator-like behavior, with a rapid but smooth crossover 
between these two phases, since a sharp MIT occurs only at lower temperatures 39 71 ). The 
critical value of the Coulomb interaction U w 5 eV is in reasonable agreement with the 
values determined spectroscopically by fitting to model calculations, and by constrained 
LDA, see 44 for details. 

To compare with the photoemission spectrum of V2O3 spectrum by Schramme et al. 14 
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Figure 11. LDA+DMFT(QMC) spectra for paramagnetic (Vo. 962 Cro. 038)203 ("ins.") and V2O3 ("met") at 
U = 4.5, 5 and 5.5 eV, and T = 0.1 eV ra 1000 K [reproduced from Ref.44]. 



and by Kim et al. 15 as well as with the X-ray absorption data by Miiller et al, the 
LDA+DMFT(QMC) spectrum of Fig. 11 is multiplied with the Fermi function at T = 0.1 
eV and Gauss-broadened by 0.05 eV to account for the experimental resolution. The the- 
oretical result for U — 5 eV is seen to be in good agreement with experiment (Fig. 12). In 
contrast to the LDA results, our results not only describe the different bandwidths above 
and below the Fermi energy (rs 6 eV and w 2 — 3 eV, respectively), but also the position 
of two (hardly distinguishable) peaks below the Fermi energy (at about -1 eV and -0.3 eV) 
as well as the pronounced two-peak structure above the Fermi energy (at about 1 eV and 
3-4 eV). While LDA also gives two peaks below and above the Fermi energy, their position 
and physical origin is quite different. Within LDA+DMFT(QMC) the peaks at -1 eV and 
3-4 eV are the incoherent Hubbard bands induced by the electronic correlations whereas in 
the LDA the peak at 2-3 eV is caused by the e a g states and that at -1 eV is the band edge 
maximum of the a\ g and states (see Fig. 10). Note that the theoretical and experimen- 
tal spectrum is highly asymmetric w.r.t the Fermi energy. This high asymmetry which is 
caused by the orbital degrees of freedom is missing in the one-band Hubbard model which 
was used by Rozenberg et al. 11 to describe the optical spectrum of V2O3. 

The comparison between theory and experiment for Cr-doped insulating V2O3 is not 
as good as for metallic V2O3, see Ref. 75. This might be, among other reasons, due to 
the different Cr-doping of experiment and theory, the difference in temperatures (which 
is important because the insulating gap of a Mott insulator is filled when increasing the 
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Figure 12. Comparison of the LDA+DMFT(QMC) spectrum 44 at U = 5eV and T = 0.1 eV xs 1000 K below 
(left Figure) and above (right Figure) the Fermi energy (at eV) with the LDA spectrum 44 and the experimental 
spectrum (left: photoemission spectrum of Schramme et al. 74 and Kim et al. ; 75 right: X-ray absorption spectrum 
ofMuller et al. 16 ). 



temperature 71 ), or the fact that every V ion has a unique neighbor in one direction, i.e., the 
LDA supercell calculation has a pair of V ions per unit cell. The latter aspect has so far not 
been included but arises naturally when one goes from the simplified calculation scheme 
described in Section 2.7 (and employed in the present Section with different self-energies 
for the dig and e g bands) to a full Hamiltonian calculation. 

Particularly interesting are the spin and the orbital degrees of freedom in V2O3. From 
our calculations, 44 we conclude that the spin state of V2O3 is S = 1 throughout the Mott- 
Hubbard transition region. This agrees with the measurements of Park et a/. 78 and also 
with the data for the high-temperature susceptibility. 79 But, it is at odds with the 5=1/2 
model by Castellani et al. so and with the results for a one-band Hubbard model which 
corresponds to S = 1/2 in the insulating phase and, contrary to our results, shows a sub- 
stantial change of the local magnetic moment at the MIT. 71 For the orbital degrees of 
freedom we find a predominant occupation of the orbitals, but with a significant ad- 
mixture of d\g orbitals. This admixture decreases at the MIT: in the metallic phase we 
determine the occupation of the (a\ g , e gl , e^ 2 ) orbitals as (0.37, 0.815, 0.815), and in the 
insulating phase as (0.28, 0.86, 0.86). This should be compared with the experimental 
results of Park et al? % From their analysis of the linear dichroism data the authors con- 
cluded that the ratio of the configurations e g e g :e g ai g is equal to 1:1 for the paramagnetic 
metallic and 3:2 for the paramagnetic insulating phase, corresponding to a one-electron 
occupation of (0.5,0.75,0.75) and (0.4,0.8,0.8), respectively. Although our results show a 
somewhat smaller value for the admixture of a\ g orbitals, the overall behavior, including 
the tendency of a decrease of the a lg admixture across the transition to the insulating state, 
are well reproduced. 

In the study above, the experimental crystal parameters of V2O3 and 
(Vo.962Cr . 038)203 have been taken from the experiment. This leaves the question 
unanswered whether a change of the lattice is the driving force behind the Mott transition, 
or whether it is the electronic Mott transition which causes a change of the lattice. For 
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another system, Ce, we will show in Section 5 that the energetic changes near a Mott 
transition are indeed sufficient to cause a first-order volume change. 

5 The Cerium volume collapse: An example for a 4 /-electron system 

Cerium exhibits a transition from the 7- to the a-phase with increasing pressure or de- 
creasing temperature. This transition is accompanied by an unusually large volume change 
of 15%, 81 much larger than the 1-2% volume change in V2O3. The 7-phase may also be 
prepared in metastable form at room temperature in which case the reverse 7-a transition 
occurs under pressure. 82 Similar volume collapse transitions are observed under pressure 
in Pr and Gd (for a recent review see Ref. 83). It is widely believed that these transitions 
arise from changes in the degree of 4/ electron correlation, as is reflected in both the Mott 
transition 84 and the Kondo volume collapse (KVC) 85 models. 

The Mott transition model envisions a change from itinerant, bonding character of the 
4/-electrons in the a-phase to non-bonding, localized character in the 7-phase, driven 
by changes in the 4/-4/ inter-site hybridization. Thus, as the ratio of the 4/ Coulomb 
interaction to the 4/-bandwidth increases, a Mott transition occurs to the 7-phase, similar 
to the Mott-Hubbard transition of the 3d-electrons in V 2 3 (Section 4). 

The Kondo volume collapse 85 scenario ascribes the collapse to a strong change in the 
energy scale associated with the screening of the local 4/ -moment by conduction electrons 
(Kondo screening), which is accompanied by the appearance of an Abrikosov-Suhl-like 
quasiparticle peak at the Fermi level. In this model the 4/-electron spectrum of Ce would 
change across the transition in a fashion very similar to the Mott scenario, i.e., a strong re- 
duction of the spectral weight at the Fermi energy should be observed in going from the ot- 
to the 7-phase. The subtle difference comes about by the 7-phase having metallic /-spectra 
with a strongly enhanced effective mass as in a heavy fermion system, in contrast to the /- 
spectra characteristic of an insulator in the case of the Mott scenario. The /-spectra in the 
Kondo picture also exhibit Hubbard side-bands not only in the 7-phase, but in the a-phase 
as well, at least close to the transition. While local-density and static mean-field theories 
correctly yield the Fermi-level peaks in the /-spectra for the a-phase, they do not exhibit 
such additional Hubbard side-bands, which is sometimes taken as characteristic of the "a- 
like" phase in the Mott scenario. 84 However, this behavior is more likely a consequence 
of the static mean-field treatment, as correlated solutions of both Hubbard and periodic 
Anderson models exhibit such residual Hubbard side-bands in the a-like regimes. 86 

Typically, the Hubbard model and the periodic Anderson model are considered as 
paradigms for the Mott and KVC model, respectively. Although both models describe 
completely different physical situations it was shown recently that one can observe a sur- 
prisingly similar behavior at finite temperatures: the evolution of the spectrum and the 
local magnetic moment with increasing Coulomb interaction show very similar features as 
well as, in the case of a periodic Anderson model with nearest neighbor hybridization, the 
phase diagram and the charge compressibility. 86 ' 87 From this point of view the distinc- 
tion between the two scenarios appears to be somewhat artificial, at least at temperatures 
relevant for the description of the a-7 transition. 

For a realistic calculation of the Cerium a-7 transition, we employ the full Hamiltonian 
calculation described in Sections 2.2, 2.3, and 2.4 where the one-particle Hamiltonian 
was calculated by LDA and the 4/ Coulomb interaction U along with the associated 4/ 
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Figure 13. Evolution of the 4/ spectral function A{ui) with volume at T = 0.136 eV (u> = corresponds to 
the chemical potential; curves are offset as indicated; At = O.lleV -1 ). Coinciding with the sharp anomaly 
in the correlation energy (Fig. 14), the central quasiparticle resonance disappears, at least at finite temperatures 
[reproduced from Ref. 51]. 



site energy shift by a constrained LDA calculation (for details of the the two independent 
calculations presented in the current Section see Refs. 51,83 and Ref. 50). We have not 
included the spin-orbit interaction which has a rather small impact on LDA results for Ce, 
nor the intra-atomic exchange interaction which is less relevant for Ce as occupations with 
more than one 4/-electron on the same site are rare. Furthermore, the 6s-, 6p-, and 5d- 
orbitals are assumed to be non-interacting in the formalism of Eq. (13), Section 2.3. Note, 
that the 4/ orbitals are even better localized than the 3d orbitals and, thus, uncertainties 
in U are relatively small and would only translate into a possible volume shift for the 
a-7-transition. 

The LDA+DMFT(QMC) spectral evolution of the Ce 4/-electrons is presented in 
Fig. 13. It shows similarities to V2O3 (Fig. 11, Section 4): At a volume per atom 
V = 20 A 3 , Fig. 13 shows that almost the entire spectral weight lies in a large quasi- 
particle peak with a center of gravity slightly above the chemical potential. This is similar 
to the LDA solution, however, a weak upper Hubbard band is also present even at this 
small volume. At the volumes 29 A 3 and 34 A 3 which approximately bracket the a-7 tran- 
sition, the spectrum has a three peak structure. Finally, by V = 46 A 3 , the central peak 
has disappeared leaving only the lower and upper Hubbard bands. However, an important 
difference to V2O3 is that the spd-spectrum shows metallic behavior and, thus, Cerium 
remains a metal throughout this transition monitored by a vanishing 4/ quasiparticle reso- 
nance. 

To study the energetic changes associated with the rapid change of the quasiparticle 
weight at the Fermi energy, we calculate the DMFT energy per site for the model Hamil- 
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tonian (13) 



£dmft = w ]T Tr(^ DA (k)G k (io;„)) 




+ U f d. 



(47) 



Here, Tr denotes the trace over the 16 x 16 matrices, T the temperature, Nk the number 
of k points, Gk the Green function matrix w.r.t. the orbital indices, H® DA (k) the LDA 
one -particle matrix Eq. (17), and 



is a generalization of the one-band double occupation for multi-band models. 

Fig. 14a shows our calculated DMFT(QMC) energies -Edmft as a function of atomic 
volume at three temperatures relative to the paramagnetic Hartree Fock (HF) energies 
-Epmhf [of the Hamiltonian (13)], i.e., the energy contribution due to electronic corre- 
lations. Similarly given are the polarized HF energies which reproduce £?dmft at large 
volumes and low temperatures. With decreasing volume, however, the DMFT energies 
bend away from the polarized HF solutions. Thus, at T — 0.054 eV ss 600 K, a region of 
negative curvature in -Edmft — £pmhf is evident within the observed two phase region 
(arrows). 

Fig. 14b presents the calculated LDA+DMFT total energy E tot (T) = -E L da(T) + 
Et>mft{T)— £'mLDA(T 1 ) where -E m LDA is the energy of an LDA-like solution of the Hamil- 
tonian (13). 88 Since both £?lda and -Epmhf — fimLDA have positive curvature through- 
out the volume range considered, it is the negative curvature of the correlation energy in 
Fig. 14a which leads to the dramatic depression of the LDA+DMFT total energies in the 
range V = 26-28 A 3 for decreasing temperature, which contrasts to the smaller changes 
near V= 34 A 3 in Fig. 14b. This trend is consistent with a double well structure emerging 
at still lower temperatures (prohibitively expensive for QMC simulations), and with it a 
first-order volume collapse. This is in reasonable agreement with the experimental volume 
collapse given our use of energies rather than free energies, the different temperatures, and 
the LDA and DMFT approximations. A similar scenario has been proposed recently for 
the 5-a transition in Pu on the basis of LDA+DMFT calculations, 48 which solves DMFT 
by an ansatz inspired by IPT and includes a modification of the DFT/LDA step to account 
for the density changes introduced by the DMFT. 49 

In a separate LDA+DMFT(NCA) calculation for Ce, we have obtained a number of 
physical quantities for both phases which may be compared to experimental values. 50 Var- 
ious static properties extracted from the calculations 50 and their counterparts from experi- 
ments are collected in Table 1 and show an overall fair to good agreement in the tendencies 
and, except for the susceptibility, the absolute values. Since the calculation of the magnetic 
susceptibility \ m Ref- 50 was based on simplifying assumptions, the absolute numbers 
cannot be expected to match experiment. However, the general tendency and especially 
the ratio between a- and 7-Ce is in good agreement with experiment. It is interesting to 
note that the experiments predict a finite Kondo screening-scale for both phases, which ac- 
tually would point toward the KVC scenario. Finally, let us compare spectral functions for 
the 4/-states calculated with the LDA+DMFT(NCA) approach to experimental data. 91 The 
photoemission spectrum for a-Ce (upper part of Fig. 15) shows a main structure between 
3 eV and 7 eV, which is attributed to 4/ 2 final state multiplets. In the calculated spectrum 
all excitations to 4/ 2 states are described by the featureless upper Hubbard band. As a 
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Figure 14. (a) Correlation energy E DMFT — B PMH p as a function of atomic volume (symbols) and polarized 
HF energy -Baphp — £pmhf (dotted lines which, at large V, approach the DMFT curves for the respective 
temperatures); arrows: observed volume collapse from the a- to the 7-phase. The correlation energy sharply 
bends away from the polarized HF energy in the region of the transition, (b) The resultant negative curvature 
leads to a growing depression of the total energy near V = 26-28 A 3 as temperature is decreased, consistent with 
an emerging double well at still lower temperatures and thus the ct-7 transition. The curves at T = 0.544 eV 
were shifted downwards in (b) by —0.5 eV to match the energy range [reproduced from Ref. 51]. 





a-Ce Theo 




7 _ Ce Thco 


7 _ Ce ny,9u 


Pa 


0.126 


0.1558 


0.0150 


0.0426 


Pi 


0.829 


0.8079 


0.9426 


0.9444 


P2 


0.044 


0.0264 


0.0423 


0.0131 


n f 


0.908 


0.8... 0.861 


1.014 


0.971... 1 


T K , [K] 


1000 


945... 2000 


30 


60 ... 95 


X, [1CT 3 emu /mol] 


1.08 


0.53... 0.70 


24 


8.0... 12 



Table 1 . Comparison between LDA+DMFT(NCA) calculated parameters for both a- and 7-phase at T = 580 K 
and experimental values 89,90 [reproduced from Ref. 50]. Pq, Pi and P2 are partial probabilities for an empty, 
singly and doubly occupied 4/-state, nf is the /-electron occupancy, Tk the estimated Kondo temperature, and 
X the magnetic susceptibility. 



consequence of the simplified interaction model all doubly occupied states are degenerate. 
This shortcoming in our calculation is responsible for the sharply peaked main structure. 
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Energy (eV) 

Figure 15. Comparison between combined photoemission" and BIS 92 experimental (circles) and theoretical 
(solid line) /-spectra for a- (upper part) and 7-Ce (lower part) at T = 580 K. The relative intensities of the 
BIS and photoemission portions are roughly for one 4/ electron. The experimental and theoretical spectra were 
normalized and the theoretical curve was broadened with resolution width of 0.4 eV. In the insets a comparison 
between RIPES 93 experimental (circles) and theoretical (solid line) /-spectra is given. The experimental and 
theoretical data were normalized and the theoretical curve was broadened with broadening coefficient of 0.1 eV 
[reproduced from Ref. 50]. 

The neglected exchange interaction would produce a multiplet structure, which would be 
closer to the experiment. The experimental peak at about 0.5 eV is attributed to two 4/ 1 
final states, which are split by spin-orbit coupling. The calculated /-spectrum shows a 
sharp quasiparticle or Kondo resonance slightly above the Fermi energy, which is the re- 
sult of the formation of a singlet state between /- and conduction states. We thus suggest 
that the spectral weight seen in the experiment is a result of this quasiparticle resonance. 
Since we did not yet include spin-orbit coupling in our model, we cannot observe the men- 
tioned splitting of the resonance. However, as it is well known, 94 the introduction of such a 
splitting would eventually split the Kondo resonance. If we used the experimentally deter- 
mined value of about 0.3 eV for the spin-orbit splitting, 92 the observed resonance of width 
0.5 eV would indeed occur in the calculations. In the lower part of Fig. 15, a comparison 
between experiment and our calculation for 7-Ce is shown. The most striking difference 
between lower and upper part of Fig. 15 is the absence of the Kondo resonance in the high 
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temperature phase (7-Ce; transition temperature 141 K &1 ) which is in agreement with our 
calculations. 

In the insets of Fig. 15, our results for the non-occupied states in the /-density are 
compared with RIPES data. 93 The calculated /-spectra were multiplied by the Fermi-step 
function and broadened with an Lorentzian of the width 0.1 eV in order to mimic the exper- 
imental resolution in the theoretical curves. Here, as above the theoretical overestimation 
of the sharpness of the upper Hubbard band is a consequence of the simplified local inter- 
action and thus of the missing multiplet structure of the 4/ 2 -final states. The main feature 
of the experimental spectra, i.e., a strong decrease of the intensity ratio for Kondo reso- 
nance and upper Hubbard band peaks from a- to 7-Ce, can also be seen in the theoretical 
curves of Fig. 15 as well as in the study presented in Fig. 13. A more thorough comparison 
of these two independent LDA+DMFT(NCA) and LDA+DMFT(QMC) studies remains to 
be done. 

6 Conclusion and Outlook 

In this paper we discussed the set-up of the computational scheme LDA+DMFT which 
merges two non-perturbative, complementary investigation techniques for many-particle 
systems in solid state physics. LDA+DMFT allows one to perform ab initio calculations 
of real materials with strongly correlated electrons. Using the band structure results cal- 
culated within local density approximation (LDA) as input, the missing electronic corre- 
lations are introduced by dynamical mean-field theory (DMFT). On a technical level this 
requires the solution of an effective self-consistent, multi-band Anderson impurity prob- 
lem by some numerical method (e.g. IPT, NCA, QMC). Comparison of the photoemission 
spectrum of Lai_ x Sr x Ti0 3 calculated by LDA+DMFT using IPT, NCA, and QMC re- 
veal that the choice of the evaluation method is of considerable importance. Indeed, only 
with the numerically exact QMC quantitatively reliable results are obtained. The results 
of the LDA+DMFT(QMC) approach were found to be in very good agreement with the 
experimental photoemission spectrum of La .94Sr .06TiO3. 

We also presented results of a LDA+DMFT(QMC) study 44 of the Mott-Hubbard metal- 
insulator transition (MIT) in the paramagnetic phase of (doped) V2O3. These results 
showed a Mott-Hubbard MIT at a reasonable value of the Coulomb interaction U w 5eV 
and are in very good agreement with the experimentally determined photoemission and 
X-ray absorption spectra for this system, i.e., above and below the Fermi energy. In 
particular, we find a spin state S = 1 in the paramagnetic phase, and an orbital admix- 
ture of e^e^ and e^aig configurations, which both agree with recent experiments. Thus, 
LDA+DMFT(QMC) provides a remarkably accurate microscopic theory of the strongly 
correlated electrons in the paramagnetic metallic phase of V2O3. 

Another material where electronic correlations are considered to be important is 
Cerium. We reviewed our recent investigations of the Ce a-j transition, based on 
LD A+DMFT(QMC) 5 1 and LDA+DMFT(NCA) 50 calculations. The spectral results and 
susceptibilities show the same tendency as seen in the experiment, namely a dramatic re- 
duction in the size of the quasiparticle peak at the Fermi level when passing from the ot- 
to the 7-phase. While we do not know at the moment whether the zero-temperature quasi- 
particle peak will completely disappear at an even larger volume (i.e., in a rather Mott-like 
fashion) or simply fade away continuously with increasing volume (i.e., in a more Kondo- 
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like fashion), an important aspect of our results is that the rapid reduction in the size of 
the peak seems to coincide with the appearance of a negative curvature in the correlation 
energy and a shallow minimum in the total energy. This suggest that the electronic corre- 
lations responsible for the reduction of the quasiparticle peak are associated with energetic 
changes that are strong enough to cause a volume collapse in the sense of the Kondo vol- 
ume collapse model, 85 or a Mott transition model 84 including electronic correlations. 

At present LDA+DMFT is the only available ab initio computational technique which 
is able to treat correlated electronic systems close to a Mott-Hubbard MIT, heavy fermions, 
and /-electron materials. The physical properties of such systems are characterized by the 
correlation-induced generation of small, Kondo-like energy scales which require the ap- 
plication of genuine many-body techniques. The appearance of Kondo-like energy scales 
in strongly correlated systems leads to several experimentally relevant consequences. One 
of the most important features is the enhancement of the quasiparticle mass m* (i.e., the 
decrease of the quasiparticle residue Z). This phenomenon can be observed as an enhance- 
ment of the coefficient 7 in the specific heat. Another important characteristic is the Wil- 
son ratio between 7 and the Pauli spin susceptibility \- Future LDA+DMFT investigations 
will determine these quantities for real systems, as well as the optical conductivity, phase- 
diagrams, the local vertex function, and various susceptibilities. 

LDA+DMFT provides, at last, a powerful tool for ab initio investigations of real mate- 
rials with strong electronic correlations. Indeed, LDA+DMFT depends on the input from 
both band structure theory and many-body approaches. Hence, for this computational 
scheme to be entirely successful in the future two strong and vital communities will finally 
have to join forces. 
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